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Abstract. Hydrogen majority plasmas will be used in the initial non-activated phase 
of ITER operation. Optimizing ion cyclotron resonance heating (ICRH) in such 
scenarios will help in achieving H-mode in these plasmas. Past JET experiments 
with the carbon wall revealed a significant impact of intrinsic impurities on the 
ICRH performance in ( 3 Hc)-H plasmas relevant for the full-field initial ITER phase. 
High plasma contamination with carbon impurities resulted in the appearance of 
a supplementary mode conversion layer and significant reduction in the transition 
concentration of 3 He minority ions, defined as the concentration at which the change 
from minority heating to mode conversion regime occurs. In view of the installation 
of the new ITER-like wall at JET, it is important to evaluate the effect of Be and 
W impurities on ICRH scenarios in ( 3 Hc)-H plasmas. In this paper, an approximate 
analytical expression for the transition concentration of 3 Hc minority ions is derived 
as a function of plasma and ICRH parameters, and accounting for typical impurity 
species at JET. The accompanying ID wave modeling supports the analytical results 
and suggests a potential experimental method to reduce 3 He level needed to achieve 
a specific heating regime by puffing a small amount of 4 He ions additionally to ( 3 He)-H 
plasma. 



PACS numbers: 52.25 Fi, 52.25 Ya, 52.55 Fa 
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1. Introduction 

Ion cyclotron resonance heating (ICRH) has been used successfully for bulk ion and 
electron heating in tokamaks, and is foreseen as one of the additional heating systems 
to be installed at ITER. At the initial stage of ITER operation predominantly hydrogen 
(H) or helium-4 ( 4 He) plasma will be used to minimize the activation of the tokamak 
components pQ. The existing scalings suggest the H-mode power threshold to be higher 
by a factor of two for hydrogen plasmas than for deuterium, and therefore the access 
to H-mode of operation is not assured for H plasmas in ITER with the heating powers 
that will be available [2]. Thus, it is particularly important to optimize the efficiency 
of radio frequency (RF) heating for the scenarios relevant for this experimental stage of 
ITER. 

It is well known that for efficient ICRH heating at the fundamental ion cyclotron 
(IC) frequency plasmas consisting of at least two different ion species need to be 
used. The concentration of one ion species (majority) is usually much higher than 
the concentration of another species (minority ions). Depending on the relative 
concentrations of the species, two heating regimes are usually identified: minority 
ion heating (MH) and mode conversion (MC). MH requires relatively low minority 
concentrations, less than some critical value [3]. In this regime the majority ions assure 
favourable polarization of the fast Alfven wave (FW) launched by the ICRH antenna 
at the region of the fundamental cyclotron resonance of minority ions, which absorb 
the RF energy and transfer it to bulk plasma ions and electrons via Coulomb collisions. 
Whether indirect bulk ion or electron heating dominates, depends on the ratio between 
the minority tail energy and the critical energy, E crit [3]. When the averaged energy of 
the fast minority ions is above E cr n electrons are predominantly heated by collisions with 
the fast ions, whereas for the opposite case indirect bulk ion heating is observed. With 
the gradual increase in the minority concentration the MH efficiency reduces and at large 
enough minority concentrations plasma heating via MC becomes dominant. This regime 
is characterized by a partial conversion of the FW to the short wavelength modes, ion 
Bernstein wave (IBW) and ion cyclotron wave (ICW), at the MC layer [5]. The converted 
wave is commonly strongly absorbed by electrons within a narrow spatial region on much 
shorter time scale than the characteristic time for indirect bulk plasma heating via MH. 
Thus, as one of the methods to identify experimentally whether MH or MC heating 
occurs, is to study electron temperature response to ICRH power modulation 0,17]. 
Beyond plasma heating itself, the MC regime has a number of promising applications 
in present-day and future fusion machines, e.g. driving non-inductive current and 
generating plasma rotation [B14T3"]. 

A number of experiments were performed at JET and ASDEX-U aimed at studying 
various ICRH heating schemes in H majority plasmas with helium-3 ( 3 He) minorities 
that can be used in the non-activated phase of ITER [T3HT8]. Both ( 3 He)-H schemes 
relevant for the half-field (Bq = 2.65 T) and full-field (Bq = 5.3 T) operation phase 
of ITER were tested recently. For the frequency range designed for the ITER ICRH 
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system (/ = 40 — 55 MHz) two heating scenarios are feasible for the half-field H phase: 
fundamental IC heating of H majority ions at / ~ 40 MHz and second harmonic heating 
of 3 He minority ions at / ~ 53 MHz [18]. However, both of these scenarios have relatively 
low single-pass absorption with dominant fast wave electron heating. For the full-field 
H phase fundamental ion cyclotron heating of 3 He ions is the only scheme available for 
plasma heating, which commonly has an as good heating efficiency as any of fundamental 
minority ICRH scenarios. This scheme has been tested at JET [15l[T6] by adopting the 
magnetic field (-Bo ~ 3.6 T) and ICRH frequency (/ ~ 37 MHz) to locate 3 He cyclotron 
layer in the plasma center as it will be in ITER. 

In JET experiments reported in Refs. [I~5l[T6] ( 3 He)-H plasma heating was studied 
at very low 3 He concentrations. Mode conversion was found to be reached at 3 He 
concentrations A[ 3 He] = m}i c /n e = 2 — 3%, which were substantially smaller than the 
values observed in ( 3 He)-D plasmas (X[ 3 He] ~ 10 — 15%). Such a difference is explained 
by the fact that ( 3 He)-H is the so-called 'inverted' ICRH scenario. A key feature of the 
inverted scenarios is that the minority ion species have a smaller charge-to-mass ratio 
than the majority species, (Z/A) mino < (Z/A) mgL j . For these heating scenarios the MC 
layer is located between the ICRH antenna on the low-field side (LFS) and the minority 
cyclotron resonance (the FW encounters the MC layer first), while for standard scenarios 
- like in ( 3 He)-D plasmas - the minority cyclotron layer is located between the MC layer 
and the LFS antenna. 

Those experiments highlighted an essential effect of impurities in the inverted 
( 3 He)-H scenario. Three important issues due to plasma dilution with the carbon 
impurities were outlined (note that the experiments [151418] were carried at JET with the 
inner vessel covered with carbon tiles). First, the heating region was found to be shifted 
appreciably away from where it was expected for pure plasma. Second, MC heating was 
complicated further through the appearance of the supplementary MC layer associated 
with carbon (C) impurities. The efficiency of ICRH heating in ( 3 He)-H JET plasmas 
with multiple MC layers was addressed in detail in Ref. [17] for the extended range 
of 3 He concentrations, showing in addition the complexity of the real time control of 
the minority level and, thus, difficulty in controlling the location of the MC layer in 
such plasmas. The third direct effect of carbon impurities on ICRH performance was 
the reduction of the transition 3 He minority concentration, X cr it[ 3 He]. Full- wave ICRH 
simulations have shown that for the plasma without the carbon the transition from MH 
to MC heating should occur at X cr ; t [ 3 He] ~ 5%, while the experimentally observed levels 
at JET were lower [T5] . 

Since August 2011 JET is operating with the new ITER-like wall, using beryllium 
(Be) and tungsten (W) as the new plasma facing materials. Therefore, it is instructive 
to assess and analyze the impact of modest amounts of first wall material impurities, 
which will enter the plasma due to the plasma-wall interaction, on the performance 
of ICRH heating. The aim of the present paper is to find a reasonable estimate for 
A cr it[ 3 He] in ( 3 He)-H plasma, which corresponds to the change of the heating regime 
from MH to MC, and evaluate the effect of typical impurities at JET on that. 
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The paper is organized as follows. Section [2] shows ID numerical results for the 
dependence of the transition 3 He concentration on the plasma and ICRH parameters in 
pure ( 3 He)-H plasma. We present two equivalent approaches for the analytical estimate 
of X crit [ 3 He] in section [31 In section [4] one of the approaches is generalized further to 
account for the impurities in the plasma. Also, the dependence of X cr i t [ 3 He] on the 
concentrations of Be and other impurities is analyzed there. Based on the results of 
previous sections, in section [5] we suggest a potential method to reduce and/or control 
X cr it [ 3 He] by using additional puffing of 4 He ions to ( 3 He)-H plasma. Finally, conclusions 
are drawn in section [61 

2. Numerical results for the transition concentration of helium-3 ions in 
pure ( 3 He)— H plasma 

To analyze the wave propagation and damping dynamics in the ( 3 He)-H plasma, the ID 
ICRH full- wave code TOMCAT [19] has been used. This code solves a 12th order wave 
equation system, that guarantees a positive definite and purely resonant absorption for 
Maxwellian populations, accounting for the radial variation of the toroidal magnetic 
field and FW parallel wavenumber. It, however, omits the finite poloidal magnetic field 
effects, and thus excludes MC of the FW to the ICW. All the results reported assume 
pure excitation of the FW from the LFS of the tokamak as the imposed boundary 
conditions. TOMCAT gives scattering coefficients (reflection, transmission, conversion 
and absorption) for a single or a double transit of the FW over the plasma depending on 
a chosen radial range of integration. In contrast to usual full-wave codes, where all the 
RF power launched into the plasma is assumed to be absorbed (i.e. yielding multi-pass 
absorption), TOMCAT calculates the single- or double-pass absorption coefficients. The 
evaluation of a single- or a double-pass absorption coefficient allows for estimating the 
heating efficiency of the studied ICRH scenario qualitatively. 

Even though TOMCAT can give only qualitative results (since many effects are 
not taken into account in a ID geometry), it is helpful for understanding the global 
trends and some of the observed characteristics of the ICRH performance. TOMCAT 
simulations were used for the analysis of the past ( 3 He)-H experiments at JET and 
helped to identify the effect of C impurities in these scenarios [151 02] • 111 a recent 
paper [20], the normalized absorption coefficients given by TOMCAT and the evaluated 
deuterium transition concentration in D-T plasma were found to be in a reasonable 
agreement with the results of more sophisticated modelling with the 2D full-wave code 
TORIC. We use therefore TOMCAT modelling to check and supplement the analytical 
estimates presented in this paper. However, it is important to note that a rigorous 
treatment of the FW propagation and MC in tokamaks should be essentially based on 
the 2D or 3D full- wave modelling [21H26]. 

Figure [H shows the single-pass (left) and double-pass (right) absorption coefficients 
in ( 3 He)-H plasma, computed with TOMCAT. We consider plasma and ICRH 
parameters typical for past JET ( 3 He)-H experiments: plasma major and minor radius 
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Helium-3 concentration, X[ 3 He] (%) Helium-3 concentration, X[ 3 He] (%) 

Figure 1. Single-pass (a) and double-pass (b) absorption coefficients vs. 3 He minority 
concentration calculated with the TOMCAT code for ( 3 He)-H plasma: B$ = 3.1 T, 
/ = 32.2 MHz, n e0 = 3.2 x 10 19 m" 3 , T = 5.0 kcV, n tor = 27. 

Ro = 2.96 m and a = 0.9 m, RF frequency / = 32.2 MHz, central magnetic field 
Bq = 3.1 T, n e (r) = (n e0 — n e i)(l — (r/a) 2 ) + n e i, central density n e0 = 3.2 x 10 19 m~ 3 , 
edge density n el = 0.1 n e0 , T(r) = (T — 71) (1 — (r/a) 2 ) 1 ' 5 + T±, central electron and ion 
temperature T = 5.0 keV, edge temperature T\ = 0.1 To, FW toroidal mode number 
is taken to be n tor = 27 (k\\ = n tor /i?), typical for dipole phasing of the A2 ICRH 
antenna at JET. This choice of / and B places 3 He cyclotron resonance almost centrally, 
i?3Hc = 2.9 m. The results we present in this section are computed assuming pure plasma 
without impurities. 

For the considered parameters single-pass absorption by minority species reaches 
its maximum pi = 27.5% at X[ 3 He] = 4.2%. Electron absorption at this 
minority concentration is only p e = 2.7%. Ion absorption starts to degrade with 
increasing X[ 3 He], whereas electron heating via mode conversion, in contrast, increases. 
At X[ 3 He] = 5.9% minority ion heating is balanced by electron heating, and we refer to 
that X[ 3 He] as a transition concentration, and denote it with X cr i t [ 3 He]. The maximum 
electron heating p e = 25.6% is reached at A[ 3 He] = 8.1%, and starts to decrease for 
higher X[ 3 He] in agreement with the Budden theory for the isolated MC layer [27] . 

Figure d^b) shows the absorption power fractions transferred in a double sweep 
of the FW in the plasma. The FW is followed from the incidence side (LFS) and is 
allowed to reflect once on the high-field side (HFS) cutoff. It exhibits oscillations in the 
double-pass absorption by electrons and significant increase in the absorption by 3 He 
ions in the MH regime. These are due to the multiple FW reflections in the plasma 
and constructive/destructive interference undergone by the reflected waves [2TJ |[28T[3Tj] . 
Electron heating via mode conversion starts to dominate over minority ion heating 
at X[ 3 He] = 6.4%. The total phase difference, which defines the resulting double- 
pass absorption by ions and electrons, includes the terms due to the FW reflection 
on the minority cyclotron resonance and MC layer, respectively. In general, it is 
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different for minority and MC heating, and thus the discussed transition concentration 
X crit [ 3 He] differs somewhat from the value calculated for a single transit of the FW in the 
plasma. Since the FW interference complicates the analysis of the wave propagation and 
absorption, and the main emphasis of this paper is to sort out the effect of impurities on 
X crit [ 3 He], the numerical results presented in the paper rely on treating the propagation 
of the single incident wave and ignoring possible additional FW reflection from the 
HFS cutoff and supplementary MC layers (single-pass TOMCAT calculations) unless 
otherwise stated. 

ICRH system on JET covers a frequency range from 23 to 57 MHz, providing access 
to a large number of scenarios for a wide range of magnetic fields [31]. In ( 3 He)-H ICRH 
experiments reported in [T5l[T7] different RF frequencies were used for plasma heating, 
/ ~ 37 MHz and / ~ 32 MHz, respectively. Figure [2(a) presents X cr it[ 3 He] as a function 
of RF generator frequency /, but for the fixed //-Bo ratio to keep the same location of 
the minority cyclotron resonance in the plasma, -R3 He . The transition concentration 
XcntpHe] decreases if operating at higher RF frequency and magnetic field: while 
for baseline case shown in Fig. 1 (f/B = 32.2/3.1) transition from MH to MC was 
reached at X[ 3 He] = 5.9%, X crit [ 3 He] decreases to 5.1% if choosing f/B = 37.4/3.6. 
X cr i t [ 3 He] follows approximately a l//-dependence, which is clearly seen by comparing 
the numerical results with the fitting curve. It suggests that the transition concentration 
is connected to the Doppler width of the minority IC resonance, AR oc y/2k\\v t h,mmo/u, 
which is also inversely proportional to the RF frequency. 

This interpretation is supported by Fig. [2(b), where the dependence of X crit [ 3 He] 
on the plasma temperature is depicted (all plasma species are assumed to have equal 
temperatures). The transition concentration of 3 He ions raises with the temperature 
and scales as a square root of T . In future machines like ITER, where higher plasma 
temperatures are expected for the tokamak operation, the transition to MC in ( 3 He)-H 
plasmas is to occur at higher 3 He concentrations than in JET. Figure[2(b) also shows that 
X[ 3 He], at which single-pass minority ion absorption is maximized (triangles), increases 
with To. Due to increased 3 He demand and industrial consumption, the typical market 
price for 3 He has raised from $100 - $200 per liter to $2,000 per liter in recent years [32J. 
Along with the fact that the plasma volume in ITER is almost 10 times larger than in 
JET, this increases significantly the operational costs for using 3 He in future fusion 
devices. The present paper discusses the possibility to retune ICRH scenarios involving 
helium-3 ions to minimize 3 He concentrations needed for MH and MC heating. 

As follows from Fig. [2(c), the transition concentration increases with the FW 
toroidal wavenumber, which is consistent with the increase in the Doppler broadening 
of the cyclotron resonance for higher kn. However, the calculated values of X crit [ 3 He] 
are less sensitive to the value of n tor than the estimates yielded from the linear fitting. 
A moderate increase in X CT \ t [ 3 He] is also observed if increasing the central plasma density, 
as Fig. [2j(d) illustrates. The reason why X CTit [ 3 He] raises with n e0 will be outlined in the 
next section, where we present simplified analytical models that help us to grasp the 
main physical parameters affecting AT crit [ 3 He]. 
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Figure 2. Transition concentration of 3 He ions as a function of RF frequency 
(//-Bo fixed) (a), central plasma temperature (b), FW toroidal wavenumber (c) and 
central plasma density (d). 



3. Analytical estimates for the transition concentration of helium-3 ions in 
pure ( 3 He)— H plasma 

ICRH heating of fusion plasmas relies on the transport of the energy by the fast waves 
from the edge of the plasma to the core. The FW propagating at the equatorial plane 
of a tokamak is fairly well described by the dispersion relation 

2 (cl -wg)(e R -njj) 

where n\\ = ck\\/(2irf) is the refractive index parallel to the equilibrium magnetic field, 
and / is the RF generator (antenna) frequency. Omitting the poloidal magnetic field, 
the FW parallel wavenumber varies radially as k\\ = n tm /R\ the dominant n tor in the 
launched FW spectrum is determined by the antenna geometry and chosen antenna 
phasing. In Eq. ([T]), e$, £l and are the plasma dielectric tensor components in the 
notation of Stix ^j, which in a cold-plasma limit and for the ICRH frequency range are 
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given by 



9 2 
U)t ^ ^pi 



*L«l + -f + E ,,v ( 2 ) 



where the summation is to be taken over all ion species constituting the plasma, and 
Up S and u cs are the species' plasma and cyclotron frequencies, respectively. 

The FW resonance condition es = njj defines the location of the MC layer, often 
called also as the ion-ion hybrid (IIH) resonance. Hot plasma theory resolves this 
resonance and bends it into a confluence. At this layer the FW is converted partially 
to a short wavelength mode, the IBW and the ICW, depending on the relation between 
the plasma temperature and poloidal magnetic field. The MC layer is accompanied 
closely to the LFS by the left-hand polarized cutoff (L-cutoff ) defined by the condition 
= tv». In Appendix A of the paper we present the derivation of the formulas for 
resonant frequencies, u$ and cjl, satisfying the conditions 

es = n 2 , 

(3) 

£l = np 

which we refer to further as the IIH and L-cutoff frequencies, respectively. 

Since the toroidal magnetic field follows 1/R dependence in tokamaks, the location 
of the layer, at which the antenna frequency u matches a certain resonant frequency uj-i 
{u = 2irf = Ui), can be calculated as follows 

^ 15.25 gg(T) 

RL ^- Ro Z^ /(MHz) ' (4) 
where cj c h is the cyclotron frequency for hydrogen ions. As an example, for the baseline 
conditions chosen (/ = 32.2 MHz, B = 3.1 T) the hydrogen IC resonance is located 
at i?H = 4.34 m, 3 He resonance - at -R3 He = 2.90 m, Be resonance - at i?Be = 1.93 m. 
Vice versa, using Eq. (j4j) we can connect the radial coordinate in the plasma, Ri to the 
corresponding resonant frequency, (Jj. 

For the discussion of the transition minority concentration and in order to 
understand how two ICRF heating regimes arise in two-ion component plasmas, we 
need to account for the kinetic response of the resonant minority ions in the es tensor 
component in the denominator of Eq. (OQ) [3j. To keep the algebra simpler, we consider 
the leading order terms in finite Larmor radius expansion of es- Then, it can be written 
as follows 

es = 1 + ^ + E f|-/^— + , (5) 

U ce i=l >2 2UJ V / 2fc||% 

where Z{£) is the plasma dispersion function, v t i = (Ti/nii) 1 ^ 2 is the thermal velocity, 
and £±ij = (ui ±u> C i)/(y/2k\\v t i). Following Ref. j3], throughout the paper indices '1' and 
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'2' refer to majority (hydrogen) and minority (helium-3) ions, respectively. Majority 
species are non- resonant, so one could use asymptotic expansion of the plasma dispersion 
function for large arguments, Z(£) ~ — 1/£. For minority ions the term proportional 
to Z(^ 12 ) is non-resonant, while the term proportional to Z(£_ 12 ) represents resonant 
minority response. In such a way, we could write es tensor element as 

S ~ u 2 -u 2 cl 2u(u + u c2 ) 2u 2 V2k llVt2 {) 

Accounting for u ~ u c2 and that the minority concentration is typically much lower 
than that of majority level, the resonance condition eg = nf, can be re-written 



7( t \ Hi ^ k \\ V t2 
^(.£-12 J - — y 



Wp 2 W 



/i 2 + k\v 2 Al 



(7) 



1 — jl 2 u 2 

where we have introduced /i = Z\A 2 jZ 2 A\, and vai is the Alfven speed corresponding 
to the majority ions. The real part of Z(j~) function is shown in Fig. [3j note that the 
argument of the plasma dispersion function is proportional to the distance from the 
minority cyclotron resonance layer (R 2 = R + x 2 ) 

\^2k\\v t2 R2 

The right-hand side of Eq. (J7J) depends on the minority concentration X 2 = n 2 /n e 
via u 2 2 in the denominator, but is independent of £_i2 and thus is represented by 
a horizontal line in Fig. [3j Note that the second term in Eq. (CO) is small compared to 
the first for the typical experimental conditions. The resonance and MC occur provided 
there is an intersection of the horizontal line with the curve for Z{^) function. For the 
standard ICRH scenarios \i < 1, while for the inverted - \i > 1; therefore, the resonance 
condition for MC can be fulfilled only on the HFS and LFS of the minority IC resonance, 
respectively. 

In Fig. [3] we have depicted 4 horizontal lines corresponding to different minority 
concentrations. When X 2 is very small, there is no intersection of the horizontal line 
(dotted) with Re[Z(£)] and thus the resonance condition for MC does not occur. The 
FW is absorbed by minority ions due to the imaginary part of Z(l;_i 2 ) and gives rise to 
the MH regime. When, on the contrary, X 2 is very high, there is always an intersection 
at |£_i2| > 5, for which the imaginary part of Z(£) is negligible (dashed line). This is 
a MC regime, when the absorption is defined by the local resonance and there is almost 
no minority damping. When the minority concentration is such that the horizontal curve 
is tangential to the real part of Z(£), then MC starts to occur (dash-dotted line). Since 
the maximum absolute value of Re[Z(£)] is approximately unity, then we can derive the 
corresponding minority concentration to be [3] 



(Wesson) 



2 tu A 2 Z\ 



— /i 2 | u 2 



(9) 



The difference between Eq. ([9]) and the result presented in Ref. [3] is '±' sign for the 
second term: plus/minus sign should be taken for the standard/inverted ICRH scenario, 
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Figure 3. MC occurs when there is an intersection of the horizontal line, which 
represents the right-hand side of Eq. ([7]) , with the real part of the plasma dispersion 
function, Z(£_i2). For the standard/inverted ICRH scenarios the resonance condition 
es = riu is fulfilled to the HFS/LFS of the minority IC layer, respectively. 

respectively. The physical interpretation for this fact is given in the Appendix B 
of the paper. Since v\ x oc l/n e , for the inverted ( 3 He)-H scenario the transition 
minority concentration slightly increases with the plasma density as it was shown in 
Fig. E(d). Note that criterion expressed by Eq. ([9]) was first indicated by Takahashi [33] . 
In Ref. [3l], Lashmore-Davies et al generalized this expression for the cases of degenerate 
resonances in (H)-D plasmas (the majority second harmonics coincides with the minority 
fundamental) and a single-species second harmonic heating. 

We can obtain the same result as that given by Eq. (Q using the different approach. 
Apart from giving additional physical insight, this approach will make it easier to include 
the effect of multiple impurities in the model. Qualitatively the transition from the 
minority heating to mode conversion heating can be explained as follows (see Fig. H]). 
Minority cyclotron resonance has a finite Doppler width, AR = p y/2k^v t 2R/co, where 
the numerical coefficient po is of the order of unity. Let denote S as a distance between the 
cold-plasma IIH resonance and the minority IC layer. For small minority concentrations 
the IIH layer is located within the Doppler broadened IC resonance (shaded area in 
Fig. HJ), and minority heating dominates (Fig. H(a)). For large minority concentrations 
the IIH resonance is located out of the region, where the cyclotron damping by minority 
ions is important, such that 5 > AR, and electron heating via mode conversion will 
become the main absorption mechanism (Fig. IU(b)). As noted in Ref. [12j, the transition 
from MH to MC is reached when 5 = AR, i.e. when the mode conversion layer passes 
through the broadened IC resonance. 
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Figure 4. Transition from MH (a) to MC (b) heating regime occurs when the IIH 
resonance passes through the Doppler broadened minority cyclotron layer (the scheme 
corresponds to the inverted ICRH scenario, for which the IIH resonance is to the LFS 
of the minority IC layer). 



For the tangential case discussed in Ref. [3], the intersection of the curves occurs at 



12 



1. At this minority concentration, MC only starts to occur and is much weaker 
than minority ion damping. MC is expected to have the same absorption strength 
as the ion absorption for higher X 2 , when |Re[Z(£_i2)]| — 0.5 (solid line in Fig. |3]). 
The argument of the plasma dispersion function, which is formally the constant po 
(cf. Eq. (El)), for this case equals to £_i2 — ±2.3. According to Eq. (BJ, the radial 
coordinate R 2 ± Ai? corresponds to the normalized frequency (a) = u/u c -u) 



u A = Z 2 l±Po y/2h\v t2 /uj 



(10) 



where Z\ = Zi/Ai denotes the ratio of the charge number to the atomic mass for ion 
species. 

At the same time the normalized IIH frequency (satisfying the condition e§ = nV) 
as a function of the minority concentration and ku is given by (see Appendix A) 



Z, 



zl-zl {zl-zlf 

2Z 1 12 1Z\ 



a/2 



:n) 



where we introduced the convenient notation f = ZiXi for the fraction of electrons 
replaced by the ion species H\ The last term in Eq. ( ITTj) describes the shift of the 



IIH resonance due to finite n tor , a 
5 = AR (corresponding to a) a = &s) 
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Inn, u, 



pli 



X 



2,crit 



Po- 



will be fulfilled if 



A, 



\Z\ 



zl\ 



± a 



4irn e e 2 /mn. The condition 



(12) 
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which is the same as ^ Wcsson ) given in Eq. (except the multiplier po), but written 
with the new notation {aZ\ = k 2 v 2 A1 /u 2 ). For the tangential case, which corresponds 
to the beginning of MC, p ~ 1 and the results given by Eqs. (Q and (fT2j) are identical. 
On the other hand, to estimate the transition minority concentration po « 2.3 should 
be used, yielding for ( 3 He)-H plasma {Z\ = 1, Z 2 = 2/3) 

X crit [ 3 He] « 2.3 ^ k W Vt2 (1.2 - 2a/3) . (13) 

An additional argument for choosing the horizontal intersection level with Re[Z(£)] 
curve in Fig. [3] to be lower than unity is due to the good correspondence between the 
numerical results and the analytical formula ( ITBl) which is presented in Figs.|2]^c) and (d). 

In this section we considered pure two-ion majority /minority plasma excluding the 
presence of impurities. The critical minority concentration in diluted plasmas will be 
upshifted or downshifted X 2tCrit = M imp X 2tC rit, if accounting for impurities that provide 
the correction factor M imp . Thus, the choice of po constant affects the absolute value 
for the transition concentration, but it is actually not important to describe the relative 
influence of impurities on that, which is a subject of the next chapter. 

4. Transition concentration of helium-3 ions in ( 3 He)— H plasma accounting 
for impurities 

The reason why impurities affect the transition concentration of minority ions can easily 
be understood by noting that the location of the IIH resonance depends on the level of 
impurity contamination. In pure plasma it is described by Eq. (lllj) . but as impurities 
are accounted for additional terms arise. In Ref. [30], three- ion component plasmas were 
considered, and the correction term describing the effect of a single impurity species was 
obtained. In Appendix A of this paper, we generalize those formulas to take into account 
any number of impurity species. It is shown, that to the lowest order the contribution 
to u)g caused by impurities is described by the expression 

xr, ~ + ^)(Z 2 + ^l^imp) (Z\ ~ Z 2 )(Z 1 - Z imp ) 

dU S ~ 2^ -=2——= v = s J2jimp- (14) 

imp V^2 T £>imp) {^2 ~ ^imp) 

Note that in the corresponding formula, derived for a single impurity in Ref. [30], due 
to a misprint a 1/ (2Z 2 ) factor is missing in a term describing the impurity response. 
If we denote X 2iCr it to be the transition minority concentration in a pure plasma, and 
A"2,crit as that in plasmas with impurities, then the following approximate formula can 
be obtained 

« z 2 + Z k [ i^h + (z? 2 ~ f 1)2 *h *Z 2 + + 

(Zl-Z 2 f ~ , ^(Zl-Zi)(Zi + Z 1 Z iinp )(Z 1 -Z h ^- 



2Z 2 aj2 + Z^ 2ZUZ$-Z?) /2/imp 

imp J- v ^ imp; 



From Eq. (ITS]) it is easy to show that 

7 I e -[ _ST — ^ im p)(^2 + ^l^imp 

h/h ~ h (i + ^{z\ - ziyzjz^zi - z 2 mp ) /imp - 
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Table 1. Rough estimates for the impurity concentrations at JET equipped with the 
new ITER-like wall. 





Be 4 + 


C 6+ 


w 46+ 


Ni 26 + 


-^imp 


2.5 x 1CT 2 


1.5 x 10~ 3 


1.0 x 10~ 4 


1.0 x 10~ 4 


/imp 


0.10 


0.009 


0.005 


0.003 


AZ cS 


0.30 


0.05 


0.21 


0.07 



In the denominator of Eq. ( fl6l) . a(Zf — Z^)jZ\ can be neglected since accounting for 
that gives the contribution oc af mp , and a number of quadratic terms have already 
been omitted when deriving Eq. ( I14p . These quadratic terms are summarized in the 
Appendix A of the paper. Thus, the transition minority concentrations in pure and 
plasmas contaminated with impurities are related as follows 

A^2,crit ^ ^ ~~ '2imp)(-2'2 ^l-^imp) r 

A2,crit imp ^1\^2 _ Amp ) 

Since in ( 3 He)-H plasmas Z\ > 2>i > Z{ mp , the impurity contamination leads to the 
reduction of the transition concentration of helium-3 ions (Mj mp < 1): 

M imp = ^j^'j « 1 - 8X[Be] - 14.6X[C 6+ ] - 33.6X[W 28+ ]- 

-62.7X[W 46+ ] - 51.4X[Ni 26+ ] - ... 
Estimates for the concentrations of impurities typical at JET are summarized in Table [TJ 
The data has been taken from Refs. [35T437] . Similarly to the effect of C ions for the old 
JET wall, the presence of Be is expected to be the main contributing impurity to the 
reduction factor M imp for the new ITER-like wall. 




Figure 5. (a) Transition concentration of hclium-3 ions in ( 3 He)-H plasma as 
a function of X[Bc] for different / and Bq. (b) Relative change of X cr i t [ 3 He] is 
independent of / and decreases almost linearly with beryllium concentration. 



Although a number of simplified assumptions have been made to derive Eq. ( TT8]) . 
it gives a remarkably good correspondence with numerical values. Figure E](a) shows 
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the transition concentration of 3 He ions as a function of Be concentration for different 
/ and Bq (f / Bq fixed), calculated numerically with the TOMCAT solver. As discussed, 
AcritpHe] is smaller if operating at higher antenna frequencies. Regardless which 
RF frequency / is chosen, the transition concentration decreases with the concentration 
of Be impurities. As shown in Fig. [2(b), the relative change of X cr i t [ 3 He] due to the 
presence of Be does not depend on / and is in a very good correspondence with Eq. ([15]) . 
M imp = l-8A[Be]. 

The numerical factors, which appear in Eq. ( fl~8l) . have been tested for other 
impurities, and Fig. [6(a) shows an excellent agreement for C ions (relevant for the old 
JET wall). The values in Fig. [6] are calculated numerically with the TOMCAT solver. 
A rough estimate for Be concentration at JET, A[Be] = 2.5% (AZ c g = 0.3), should 
result in the reduction of the transition concentration of helium-3 ions by about 20%; 
it is smaller than the effect of carbon impurities reported previously (A[C 6+ ] = 2.5% 
corresponds to AZ cS = 0.75 and M imp = 0.63). Accounting for another impurities 
present in the plasma leads to a small further decrease of M imp , M imp « 0.75. Thus, 
for similar experimental conditions the transition from MH to MC is expected to occur 
at somewhat higher 3 He concentrations for the new JET wall in comparison with those 
reported for the old carbon wall. Along with the approximate formula ( fl3|) . Eq. (fl8|) can 
serve as a rough estimate for the transition concentration of helium-3 ions in ( 3 He)-H 
plasmas for JET-like plasmas. 
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Figure 6. Transition concentration of helium-3 ions in ( 3 Hc)-H plasma as a function 
of carbon (a) and helium-4 (b) concentrations. Figure (b) is for / = 32.2 MHz and 
B = 3.1 T. 



5. Suggestion for puffing 4 He ions to ( 3 He)— H plasma to control transition 
from MH to MC regime 

In view of the temperature dependence of A crit [ 3 He] shown in Fig. [2(b) and the 
arguments concerning the price of A[ 3 He], using smaller helium-3 concentrations for 
running ( 3 He)-H experiments is beneficial. Depending on the specific goals envisaged for 
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the ICRH system one might need to set the experimental conditions to reach either MH 
(X[ 3 He] < X crit [ 3 He]) or MC heating (X[ 3 He] > X crit [ 3 He]): operating in MH regime is 
preferred to increase the fraction of thermal ion heating in ( 3 He)-H plasmas, whereas MC 
could be potentially used for driving the current or local plasma control. Therefore, the 
development of the methods to decrease and possibly control the minority concentrations 
needed for efficient ICRH performance is of high importance. As discussed in the 
previous section of the paper, plasma dilution with impurities leads to the reduction 
of AT[ 3 He], at which the transition from minority heating to mode conversion occurs. 
This section discusses briefly the possibility of X crit [ 3 He] reduction by puffing additional 
gas to ( 3 He)-H mixture. 

The concentration of intrinsic Be impurities in a plasma is determined mostly by 
the level of plasma-wall interaction. It cannot be changed in a wide range and within 
a short time period that is a prerequisite to control the transition concentration of 
minority ions. However, additional puffing of 4 He ions to ( 3 He)-H plasma seems to 
potentially fulfil this goal. 4 He is a D-like species, and its price is much lower than the 
price for helium-3. As follows from Eq. (TIE]) , we expect that X crit [ 3 He] decreases with 
the concentration of helium-4 ions according to 

M imp « 1 - 8X[Be] - 4.9X[ 4 He]. (19) 

Figure [6]^b) shows the comparison of the results given by Eq. (fl9"|) and the values 
calculated numerically with the TOMCAT solver for the baseline parameters considered 
(/ = 32.2 MHz and B = 3.1 T). In plasmas without Be ions, addition of 5% and 10% 
of helium-4 ions to the ( 3 He)-H plasma decreases the transition concentration X cr i t [ 3 He] 
by 22% and 44%, respectively. Accounting for the presence of Be ions in the plasma, 
even smaller fractions of helium-4 can be used for minority ion control. As Fig. [7] 
illustrates, for plasmas including Be impurities at the level of 2%, by puffing 7% of 4 He 
one obtains similar transition concentration X crit [ 3 He] as those in pure ( 3 He)-H plasmas 
with X[ 4 He] = 10% puffing. 

If a fraction of fuel hydrogen ions are replaced by 4 He ions the effective charge Z e g 
will raise. But at the same time the level of 3 He, which also contributes to Z e g- , decreases, 
so the resulting increase of the effective charge is reasonably small, AZ c s « 0.1 — 0.15. 
Therefore it seems to be worth exploring this option in detail in separate studies. 

A similar potential method to control X cr i t [ 3 He] could be applied if operating with 
deuterium majority plasmas, in which helium-3 is partially burn down via the fusion 
reaction D + 3 He -»■ 4 He(3.6 MeV) + p(14.7MeV). By adding X[E] = 15% to ( 3 He)-D 
plasmas, the transition concentration of helium-3 ions in such plasmas can be reduced 
by a factor of ~ 0.75 — 0.8. 

6. Conclusions 

ITER will start its operation using the hydrogen majority plasmas to minimize the 
activation of the tokamak components. Since the access to H-mode in that plasmas will 
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Figure 7. Transition concentration of helium-3 ions in ( 3 He)-H plasma increases 
with raising the plasma temperature. By puffing 4 He ions, X cr j t [ 3 He] can potentially 
be reduced and/or controlled. 



not be assured with the available heating powers, it is of high importance to maximize 
the heating performance for the non-activated phase of ITER operation. ICRH is one 
of the heating systems to be used in ITER, and a number of ICRH experiments were 
performed at JET to develop and optimize heating scenarios relevant for such plasmas. 
One of the promising ICRH schemes relies on the use of the resonant minority 3 He ions 
to absorb the RF energy. 

Past JET experiments with the old wall consisting of carbon tiles highlighted 
a number of specific features when applying ICRH in ( 3 He)-H plasmas. Particularly, 
a significant effect of intrinsic carbon impurities was revealed in that experiments. The 
heating region was found to be shifted radially comparing to the heating maximum 
expected in a pure plasma, and that was attributed to the presence of impurities. Under 
certain conditions C impurities also produced a supplementary MC layer, and multiple 
MC dynamics was observed in ( 3 He)-H plasmas. In addition, the concentration of 
minority 3 He ions, at which the transition from MH to MC occurs, was reported to 
be lower than that predicted by numerical simulations if the effect of impurities was 
neglected. 

In the present paper we discuss how the transition concentration of 3 He ions depends 
on the plasma and ICRH parameters. It is shown that X crit [ 3 He] is related to the Doppler 
broadening of the minority IC layer, and thus is inversely proportional to the antenna 
frequency and increases with the plasma temperature and FW toroidal wavenumber. 
Using two equivalent analytical approaches, we generalize the formula for the transition 
minority concentration given in Ref. [3] for the case of inverted ICRH scenarios, which 
( 3 He)-H scheme belongs to. This allows us to explain a small increase of X crit [ 3 He] with 
the plasma density in ( 3 He)-H plasmas. 
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Accounting for multiple impurity species always present in the plasma, we show 
that X crit [ 3 He] decreases and scales almost linearly with the impurity concentrations. 
An analytical estimate for the relative change of X crit [ 3 He] due to impurities is derived, 
and is shown to be in a good correspondence with the numerical results. We demonstrate 
that Be is to be the main impurity species affecting A cr i t [ 3 He] in ( 3 He)-H plasmas for 
JET equipped with the new ITER-like wall. A reduction of A crit [ 3 He] by ~ 20 - 25% 
is predicted if considering typical Be and another impurity concentrations at JET. 
A possible method to reduce and control 3 He level, at which the transition from MH to 
MC is to occur, is suggested: the method relies on the additional puffing of 4 He ions to 
( 3 He)-H plasmas. We show that for X[ 4 He] = 5 — 10% the reduction of the transition 
concentration of helium-3 ions by a factor of ~ 0.6 can be expected. 

Since a number of simplifying assumptions have been made in this paper, the results 
can serve only as a qualitative estimates rather quantitative predictions. However, our 
results support the earlier experimental findings and allows to understand the basics 
of the underlined ICRH physics and interpolate the discussed effect accounting for new 
impurity species present in a plasma. Further numerical modelling and experimental 
support is needed to confirm the potential of using 4 He gas puffing for X crit [ 3 He] control 
in ( 3 He)-H plasmas. 
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Appendix A. IIH and L-cutoff frequencies accounting for multiple impurity 
species and n tor ^ 



In the IC frequency range the cold plasma tensor components are given by 

(a.i; 



, .2 2 

s ~ + 7J~ ~ ^ 77* — TJ' 



OJ 2 (J 2 - 
pi 



The sum in (1A.1I) is to be taken over all ion species, including impurities. The conditions 
es = nf\ and 6l = n 2 defining the location of the IIH resonance and L-cutoff can be re- 
written in a simpler form 

E J lZl r,2 = a > E 7 f \~ ; = a ' ( A - 2 ) 

by introducing the following notations for the ion species t V: X{ = ni/n e , f = ZiXi 
- the fraction of the replaced electrons, Zi = Zi/Ai - ratio of the charge number to 
the atomic mass; o; c h = ei?/(m H c) - the cyclotron frequency of hydrogen ions that is 
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n 



Ann e e 2 / ttih - the reference 
— 1 — u 2 /u 2 e - effective n\ 



used for the frequency normalization (oj = u/u c n), u p u 
hydrogen plasma frequency (note n e instead of n H ), n 2 = 
if accounting for the small vacuum and electron contributions to the tensor component 
es- Then, the small parameter appearing on the right-hand side of Eqs. ( 1A.2I) . which 
describes the effect of finite n tor on the location of the IIH resonance and L-cutoff, is 
given by a = (w^/wJ H ) 
By substituting fi 

to the following ones 



1 — f 2 — /imp, it can be shown that Eqs. (IA.2I) are equivalent 



imp 



{Z x - Z 2 ){ujI + Z x Z 2 



~ 2 
UJ5 



Zl 



-h 



7s 



z\) - E 



[Zi - z 



imp 



7s, 



z x z 



(A.3) 



impy 



r,2 



imp 



? 2 
■"imp 



fi 



imp 



and 



-1 + - 



Zo 



7l 



- Z 2 
a(u) h - Z\ 



h = 7l, 



El 1 



imp 



7. 

Jimp- 

•^-"imp 



(A.4) 



The right-hand side of Eqs. (IA.3[) and (IA.4[) includes contributions due to finite n tor and 
impurities. Neglecting those terms the well-known expressions for the IIH and L-cutoff 
frequencies are obtained 



^so — Z^ 



Z- 



1 - (1 - Z x jZ 2 )h 



- Z 2 /Z 1 )f 2 

z 2 2 ) . (z 1 

-J2 



z 2 )(z 1 



Z 2 ) 2 (Z 1 



8Z 2 Z 2 



3Z 2 ) 2 
J 2 i 



(A.5) 



2Z X 

^lo = Z 2 + (Zi — Z 2 )f 2 . 

For a plasma consisting of A" ion species with different Z/A ratios and if n tor ^ 0, 
there are A" different solutions for u)g and u)l satisfying Eqs. flA.31) and (IA.4[) . Thus, 
in the general case the solutions has to be evaluated numerically. There is an artificial 
solution, which appears due to a ^ (cDl — —1/a, u)| ~ —Z 1 /a), while the other 
N — 1 solutions represent IIH resonances and L-cutoffs associated with the minority 
ions and N — 2 impurity species. We are interested in the solutions for Cos and cul, which 
correspond to the minority ions and generalize Eqs. (IA.5I) . aiming for identifying the 
influence of nt or and impurities on that. 

Omitting the comprehensive algebra, we provide the formula for L-cutoff valid up 
to quadratic terms in a and f mp , including the cross-terms: 



k u af 2 (l - f 2 ) + ^12/2/5 



imp 



imp 



1 + k 21 a + k 2 2f 



imp 



imp 



(A.6) 
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where the parameters k^ are given by 
k u = (Zt - Z 2 )\ 

k\2 — {Z\ ~ Z 2 ){Zi — Z imp )/(Cj^o ~~ ^imp) 

k 21 = (Z 1 -Z 2 )(l-2f 2 ), 

[Z\ — Z imp )(Z 2 — Z- 1Ta p) I (£>lo — Z\ 



(A.7) 



^22 



'imp / 



The expression for the IIH resonance is more complicated since Eq. (1A.3|) is written for 
cDg rather than In that case, it can be shown that 



u 2 



so 



k u af 2 (l- f 2 ) + J2ki2f2f 



imp 



imp 



1 + k 21 a + J2 k ^fi 



imp 



imp 



(A. 



is a high-accuracy quadratic approximation for cug. The expansion coefficients k^ for cjg 
are given by 

fcn ^i = Z\Z 2 {Z\ — Z 2 ) 2 , 



k\ 2 Z\ — Z 2 {Z\ — Z 2 )(Zi — ^i mp )(a)g + ZiZi mp )/ (ajg — Z? 



imp / ' 



^21^1 

-z u 



Z^)(Z 1 -(Z 1 + Z 2 )f 2 



(A.9) 



-^imp)/(^S0 



imp/ 



(a) so Z imi 



X^so 



2x3 



imp / 



'^SO 



(Zi- . 

yp\Z\ + ^imp)(^-'s 

where Z x = Z x - (Z 1 - Z 2 )f 2 . Note that Eqs. flA~MA9|) are valid for an 
arbitrary number of impurity species, and include the impurity and n to r contributions 
simultaneously. 

Figure lAlT a) shows the location of the IIH resonance and L-cutoff in pure ( 3 He)-H 
plasma as a function of the FW toroidal wavenumber (the other parameters correspond 
to those used in Fig. [1]). Considering ku ^ 0, both the IIH resonance and L-cutoff 
shift towards the LFS, and the thickness of the evanescence layer gradually decreases. 
As Fig. lAlf a) clearly illustrates, this shift is almost linear in a and thus quadratic in 
n tor . Analytical approximations given by Eqs. (1A.6HA.9j) are shown in Fig. IA1I with 
dots and are in excellent agreement with the numerical results for the whole range of 
experimentally relevant n toI . The formulas (IA.6HA.9"]) are in an almost perfect agreement 
with numerical results if also impurities are accounted for. Figure IAl( b) is calculated 
for n tor = 27 and illustrates further (almost linear) shift of the MC layer and L-cutoff 
towards the LFS if the beryllium concentration is increased in a ( 3 He)-H plasma. 
In contrast to the effect of n tor , the width of the evanescence layer increases gradually 
with X[Be\. 
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Figure Al. Location of the IIH resonance and L-cutoff in ( 3 He)-H plasma 
(X[ 3 Hc] = 5.9%) as a function of the FW toroidal wavenumber (a), and beryllium 
concentration (n tor = 27) (b). 
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To the lowest order (linear in a and /imp) we have 
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(A.10) 



(AH) 



The first terms in Eqs. ( lA.llj) show that both the IIH resonance and L-cutoff are 
located close to the minority IC layer. The second term illustrates a successive shift of 
the MC and L-cutoff layers towards the HFS or the LFS (depending on Z\j Z 2 ratio) 
with increasing the concentration of minority ions. The third term corresponds to the 
contribution due to finite n tor - If is always positive and, thus, accounting for k\\ 7^ 
results in a small shift of the layers towards the LFS. Finally, the fourth term represents 
the effect of impurities on the location of the layers in a plasma; its sign, which defines 
the direction of the additional shift due to impurities, depends on the sign of the ratio 
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(Zi — Z 2 )(2i — Z imp )/(Z 2 — Z imp ). Equation ( lA.lip for oj s is used in the main text of 
the paper to evaluate the transition concentration of helium-3 ions in ( 3 He)-H plasma 
and the effect of impurities on that. 

Appendix B. On '±' sign in formulas for the transition minority 
concentration 



There is a simple physical explanation for the appearance of '±' sign in Eqs. ((9]) and 
(I12p for the transition minority concentration. As mentioned in section [3J '+' or '— ' sign 
is to be taken for the standard and inverted ICRH scenarios, respectively. In plasmas 
without impurities the IIH frequency is given by 

^1 ~ ^2 t i (^1 



Z^ 



h 



Z\f 



01 f 2 



(B.l) 



2Z 1 J * ' 2Zl 

where a oc (n 2 or /n e )(B/ f) 2 . The sign of the second term in Eq. (IB.lj) is different 
for the standard {Z\ < Z 2 ) and inverted (Z% > Z 2 ) ICRH scenarios. While for the 
standard scenarios the MC layer shifts towards the HFS with increasing the minority 
concentration, for the inverted scenarios it moves in the opposite (LFS) direction 
(Figure IB1I) . At the same time the shift of the MC layer due to finite fcy (described 
by the third term in Eq. (1B.1I) ) is always towards the LFS regardless the relation of Z/A 
ratio for majority and minority ions. Thus, for the inverted ICRH scenarios the shift of 
the MC layer due to X 2 and fey is in the same direction, and therefore smaller minority 
concentration is needed for the MC layer to pass through the Doppler broadened 
minority IC region. Vice versa, for the standard scenarios the fcy correction counteracts 
the HFS shift of the MC layer with increasing X 2 . This results in higher minority 
concentrations needed to pass the border, which marks the transition from MH to MC. 



Standard scenario, e.g. ( 3 He)-D: '+' 



(HFS) 
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— ►*nt 
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Figure Bl. The shift of the IIH resonance with increasing X 2 and kn is in 
the opposite/same direction for the standard/inverted ICRH scenarios, respectively. 
It results in a '±' sign that appears in the correction term for the transition minority 
concentration given by Eqs. ([9]) and (fT2]) . 
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